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We study magnetically stabilized nematic order for spin-one bosons in optical lattices. We 
show that the Zeeman field-driven quantum phase transitions between non-nematic Mott states and 
quantum spin nematic Mott states in the weak hopping limit are in the universality class of the 
ferromagnetic XXZ (S = 1/2) spin model. We further discuss these transitions as condensation 
of interacting magnons. The development of 0(2) nematic order when external fields are applied 
corresponds to condensation of magnons, which breaks a U(l) symmetry. Microscopically, this 
results from a coherent superposition of two non-nematic states at each individual site. Nematic 
order and spin wave excitations around critical points are studied and critical behaviors are obtained 
in a dilute gas approximation. We also find that spin singlet states are unstable with respect to 
quadratic Zeeman effects and Ising nematic order appears in the presence of any finite quadratic 
Zeeman coupling. All discussions are carried out for states in three dimensional bipartite lattices. 
PACS number: 03.75.Mn, 05.30.Jp, 75.10.Jm. 



I. INTRODUCTION 

When the dimensionless exchange coupling strength is 
strong enough, Mott states of spinful particles are known 
to develop certain spin order. In the opposite limit, quan- 
tum fluctuations usually restore the broken symmetry re- 
sulting in spin singlet states. This widely accepted belief 
however does not exclude, and furthermore implies cer- 
tain hidden fluctuating order in the symmetry restored 
states. The preexisting order appears dynamically at cer- 
tain time and length scales and in general is very relevant 
to the low energy physics. Particularly if an external 
magnetic field is applied, nontrivial spin (magnetic) or- 
der might be induced because of the coupling between 
magnetic excitations and external fields. In other words, 
an external field can stabilize a spin order in a parameter 
regime where order is absent in a zero field. 

The possible field-induced ordering usually results 
from condensation of magnetic excitations, or magnons. 
The field-induced quantum phase transitions between 
states of different magnetic correlations and possible 
magnetically stabilized order close to critical points can 
be investigated by examining magnon excitations in one 
of the phases involved in the phase transitions. 

A well-known example is the S = 1 antiferromagnctic 
spin-chain 1-7 . The ground state is the spin singlet AKLT 
state and all spin excitations are fully gapped by the Hal- 
dane gap. An applied Zeeman field along the positive 
z-direction, although it has zero coupling with the sin- 
glet ground state, couples to spin excitations and lowers 
the energy of excitations in the S = 1, S z = 1 branch be- 
cause of the Zeeman coupling. At a critical field, the zero 
momentum excitation becomes degenerate with the spin 
singlet ground state signifying a quantum phase transi- 
tion. 



The presence and nature of induced canted Neel order 
in this case therefore depend crucially on interactions 
between magnetic excitations in the AKLT phase. One 
can easily visualize that the transition is of first order if 
the interactions of condensed magnons are attractive or 
absent. Naturally, the magnetization per lattice site in 
this case jumps by a finite value at the critical point as a 
result of condensation; furthermore the canted Neel or- 
der would not appear in this case. The phase transition 
would be simply between a fully polarized state and a 
spin singlet state. 

However, if interactions between magnons are re- 
pulsive, condensation takes place continuously because 
of finite chemical potentials for repulsively interacting 
magnons. Thus the magnetization per site, which is 
proportional to the density of magnons varies continu- 
ously across the critical point and the transition is of 
second order. In this case, the resultant state has canted 
Neel order. For the S — 1 antiferromagnetic spin chain, 
numerical results show that the magnetization indeed 
varies continuously and imply that magnons have repul- 
sive interactions 7 . In fact, at higher magnetic fields, ex- 
ternal fields do induce canted Neel order in S = 1 spin 
chains (only quasi-long-range order prevails in chains). 
Condensation of magnons has also been recently studied 
in three-dimensional frustrated magnets (see for instance 
Rcf. ( 8, 9)). 

Therefore, to investigate magnetically stabilized or- 
der, it is important to understand interactions between 
magnons or magnetic excitations. Generally, microscopic 
calculations of magnon interactions are not only very dif- 
ficult but also practically impossible for low-spin systems 
because of uncontrollable approximations involved. How- 
ever, in the case we are going to examine we do evalu- 
ate the interactions microscopically in various situations; 
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therefore we believe the results about magnetically sta- 
bilized quantum spin nematic order and quantum phase 
transitions are precise in this sense. We also want to 
emphasize that magnon condensation which interests us 
in this article occurs in both high dimensional and one- 
dimensional non-frustrated optical lattices; and the one- 
dimensional limit will be treated in a separated paper. 

The purpose of this article is to understand magneti- 
cally stabilized nematic order of spin-one bosons in opti- 
cal lattices. As emphasized above, our starting point will 
be a series of Mott states with no nematic order, some of 
which also have zero coupling with external (linear) Zee- 
man fields; our main subject is to investigate the develop- 
ment of nematic order or spontaneous symmetry break- 
ing in the xy-plane when external magnetic fields are ap- 
plied along the ^-direction. In particular, we will focus 
on the nematic order close to critical points where our 
results are actually exact. Furthermore, the approach we 
employ here is believed to yield exact phase boundaries 
between nematic states and non-nematic Mott states, a 
rather remarkable conclusion thanks to a powerful map- 
ping developed below. 

Spin-correlated Mott states for spin-one bosons have 
recently attracted considerable interest 10-14 . Theoreti- 
cal works indicate that spin correlations in Mott states 
depend on three dimensionless parameters 11 ' 12 . The first 
one is the dimensionless exchange coupling r](= J ex /E s ) 
which is defined as the ratio between the exchange inter- 
cation J ex and the bare spin gap for an individual site 
E s . The second parameter is the parity P (even or odd) 
of the number of particles per site in Mott states. The 
last parameter D is the dimensionality of optical lattices 
(we assume all lattices are bipartite). 

For D = 2,3, it was argued that Mott states for all 
odd P arc ncmatically ordered. In fact irrespective of the 
exact number of particles per site, the effective Hamilto- 
nian in the small hopping limit (rj <C 1) is equivalent to 
the bilinear-biquadratic model for S = 1 spin chains 14 . 
For even P on the other hand, nematic and spin singlet 
Mott states are present for large r\ and small rj limits 
respectively. The one-dimensional physics (D = 1) is 
dominated by quantum fluctuations. Both dimcrized va- 
lence bond crystals and non-degenerate spin singlet states 
have been found. As a result of symmetry restoring, the 
low energy dynamics in Mott states is mapped into the 
even- and odd-class quantum dimer models 11 . Further- 
more, atoms have a tendency to be fractionalized into 
solitonic elementary excitations in this limit. Supcrfluid 
phases have distinct topological properties and remain to 
be fully understood. 

Responses of correlated states of spin-one bosons to 
external fields are fascinating. For nematic condensates, 
the responses are continuous. The linear coupling be- 
tween condensates and weak external fields pins the easy 
axis in the xy-plane perpendicular to external fields, and 
0(3) nematic condensates become 0(2) ones, or canted 
nematic states; quadratic coupling however, pins the ne- 
matic easy axis along the direction of coupling 19 . For 



spin singlet condensates, the magnetization jumps dis- 
continuously result of a series of level crossings be- 
tween states \S,S Z = S) and \S + 2, S z = S + 2) 15 (also 
see general discussions about condensates in Ref.( 16, 
17)). Responses of Mott states to external Zeeman fields 
and various transitions have been recently studied in a 
mean-field approach. Mott states can cither respond 
to external fields continuously like nematic condensates 
or develop interesting magnetization plateaus similar to 
charge quantization in a Mott state 18 . 

In this article, we study magnetically stabilized ne- 
matic order in optical lattices. Particularly we demon- 
strate the development of nematic order as repulsively in- 
teracting magnons in non-nematic Mott states condense. 
We investigate the induced nematic order associated with 
the spontaneous breaking of 0(2) nematic symmetry the 
magnetization and the spin wave velocity. 

In section II, we review the properties of spin singlet 
Mott states and introduce a projected nematic order pa- 
rameter for discussions on spin partially polarized states. 
For a given lattice site with two particles, we show that 
nematic order can be established if a spin singlet state is 
a superposition with a higher spin state \S = 2, S z = 2). 
The relative phase between these two states in the super- 
position determines the easy axis of the nematic order 
parameter, or the orientation of spin nematic states. In 
section III, we study the general characterization of ne- 
matic order in spin polarized Mott states; we propose a 
projected nematic order parameter which projects away 
trivial contributions from spin polarization. 

In section IV, we truncate the Hilbert space close 
to critical magnetic fields and show that the resultant 
Hamiltonian is an XXZ (S = 1/2) pscudo spin model 
in an effective field along the z-direction. We carry out 
microscopic calculations of all parameters in the effective 
ferromagnetic XXZ model. These calculations are done 
for two particles per site, four particles per site; in the 
large ./V (even) limit; the quantum rotor model studied 
in previous works is employed to facilitate calculations. 

In section V, close to critical lines and a tri-critical 
point we further study the properties of various phases of 
the XXZ model in both semi-classical approximation and 
dilute gas approximation based on the Holstein-Primakov 
boson representation. We analyze the interactions be- 
tween Holstein-Primakov bosons or magnons. We obtain 
the exact phase boundaries for the ferromagnetic XXZ 
model by investigating the instability lines of magnon 
excitations. We also discuss the relation between con- 
densation of magnons close to critical lines, the variation 
of magnetization across critical points, and the appear- 
ance of ferromagnetic order. 

In section VI, we investigate, in details, the develop- 
ment of magnetically stabilized nematic order by exam- 
ining results following the mapping to the XXZ model 
and to the Holstein-Primakov condensation problem. 

We notice that results about phase boundaries, ne- 
matic order and spin wave velocities in the critical regime 
can be obtained in a dilute gas approximation and there- 
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fore, remarkably, are exact. In section VII, we further 
study the effect of quadratic Zeeman coupling. Finally, 
in section VIII, we conclude our investigation on this sub- 
ject. 



In SSM states, (Q a j3{k)) — as a result of the rota- 
tional invariance of the wave function and thus there is 
no nematic order. The hidden fluctuating nematic or- 
der can be studied by examining higher moments. For 
N = 2, one can easily obtain the following results 



II. SPIN SINGLET MOTT STATES, 
FLUCTUATING NEMATIC ORDER AND 
PROJECTED NEMATIC ORDER 

The Hamiltonian for spin-one bosons with antiferro- 
magnetic interactions in optical lattices in an external 
field can be conveniently expressed as 10-13 ' 18 , 

H = E s ]T S I + E c J2 Pk - Mo J2 & - E § kz H z 

h k k k 

-*£>L^a + fc.c). (i) 

(ki) 

Here Vfe a (V'fea) 5 a = x,y,z are creation (annihilation) 
operators for spin-one particles in three different states 

at site k. S ka = -ie a p y ipl a ipkp and Pk = Wka^ka are 
the spin and number operators defined for each lattice 
site. It is easy to verify that [Ska, Sk'fj] — i6kk> tap-ySk-y- 
E c ^ s are "bare" charge and spin gaps studied in Rcf.( 11); 
y«o is the chemical potential. Finally, the sum over (kl) 
represents the sum over all neigboring sites. 

In Mott states, the Hamiltonian can be reduced to the 
following effective one, 

H = E s Sfc - 22 ^ kz ^ z 

k k 

-J ex ^{Qati{k)Qp a {l) + h.c); 

{kl) 



(2) 
t 2 /2E c 



Eq.(2) is valid when (<£ c and E s <C E c ; J e 
is the exchange interaction. 

The effective Hamiltonian with zero Zeeman coupling 
was obtained in a few previous works; solutions to this 
Hamiltonian have been studied in various limits. When 
the external field is absent and r\ — J ex /E s is much less 
than unity, the ground state is a spin singlet for an even 
number of particles per site (N). The spin singlet Mott 
(SSM) ground state in this limit is the product of spin 
singlets at each individual site (up to a normalization 
factor) 11 - 13 , 



yr, TT (V'fea 



1>ka 



)N/1 



#+!)! 



\vac > 



(3) 



Meanwhile, the spin fully polarized (SFP) ground state 
is the product of on-site maximally polarized states , 



SFP 



=n 



vac > 



(4) 



2 2 

g^fcfc' (Sa'pSp' a + 8 a > a 5pp> - -6 a ip>5 a p), (5) 



which indicate on-site fluctuating nematic order. More 
explicitly, one finds the amplitude of fluctuations of ne- 
matic tensor matrix elements 



(6) 



To investigate the responses of spin singlet Mott states 
or other non-nematic states, which exhibit certain fluc- 
tuating nematic order, to external fields, it is important 
to understand how nematic order can be induced by ex- 
ternal Zeeman fields. For this purpose, we focus on the 
simplest situation where J ex is zero and treat each site in- 
dependently. We would like to demonstrate the following 
important statement: nematic order appears whenever a 
spin singlet state and a polarized state are in a linear 
superposition. 

We first consider two particles at one lattice site. The 
Hilbert space is spanned by five-fold degenerate S = 2 
states and a spin singlet state. When an external field 
is applied along the z-direction, the five-fold degeneracy 
is completely lifted while the maximally polarized state 
\S = 2, S z = 2) approaches the spin singlet ground state. 
When the level crossing takes place, the spin projection 
along ^-direction jumps by 2K. It is obvious that no ne- 
matic order is induced in this simple limit and there are 
no transitions between nematic states and spin singlet 
states. 

However, at the level crossing points, one can further 
study the properties of coherent superposition of |0, 0) 
and 1 2, 2) states while the rest of states are highly excited 
ones at these crossings. Let us introduce a coherent state 
defined in the two-state subspace as 



|«) = cos - exp(-zp T) + sin - exp(*|)| |); 



|T) = \S = 2,S Z = 2) 



f 



2V2 



(V4 +^lf\vac), 



\l) = \S = 0,S x =0) = -j=1>Wa\vac). (7) 



Here the unit vector is defined as Q 
(sin 9 cos cj>, sin sin <fi, cos 9) . One can easily verify that 



(S z ) =2?icos 2 -. 



(8) 
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S z reaches the maximum when 9 = and the minimum 
when 9 = tt/2. 

Direct calculations of the usual nematic order param- 
eter Q a p defined as the expectation value of the tensor 

operator Q a fj in Eq.2 suggest that it have a nontrivial 
structure in the maximally polarized state. Namely it 
contains a) an antisymmetric tensor as a result of spin 
polarization and b) a traceless symmetric part which re- 
flects the explicit rotational symmetry breaking by the 
magnetic field but is not associated with the spontaneous 
nematic symmetry breaking in a plane perpendicular to 
the polarization. 

To discuss the nematic order in fully or partially po- 
larized states, it is therefore essential to introduce a pro- 
jected nematic order parameter Q^p(a,/3 = x,y, z), 



Qa/3 



(9) 



Note that in the projected order parameter, the com- 
ponent associated with the nematic symmetry breaking 
remains while the components associated with polariza- 
tion have been projected away. Two tensors we would 
like to project away are defined as 



Two important features in Eq.(ll) are worth emphasiz- 
ing. Firstly, the two eigenvalues correspond to ± sin 9/^/3 
and are proportional to the coherence factor in the linear 
superposition of coherent states. They are nonvanishing 
only if 9 is not zero or 7r. Therefore nontrivial nematic 
order always appears when | f) and 1 1) two states are in 
a linear superposition. 

Secondly, the eigenvector with the maximal eigen value 
represents the easy axis of nematic order. And the easy 
axis is fully characterized by the azimuthal angle of $1. 
Indeed, one finds that the easy axis in the ccy-plane is 
defined as a 2D unit vector in the xy-plane. 



(cos • 



(12) 



When varies from to 2tt, the easy axis ui rotates by w 
angle in the xy-plane. And the nematic order parameter 
is indeed a tensor constructed out of the 2D unit vector 



2 1 
Qaf3 V = -=sin9(uj a ujp - -S aP ). 



(13) 




(10) 



One finds that the nematic order only appears in the 
xy-plane perpendicular to the external fields. It is in- 
deed straightforward to show that the truncation of Q^r 
in the rry-plane (i.e. elements with a 
coherent state defined above is 



x, y only ) for a 



n Px v 



1 
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cos <p sin (f> 
sin <f) — cos ( 
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FIG. 1. Coherent states |fi) at the Bloch sphere of pseudo 
spins. All states except the north pole (9 = 0) and south pole 
(6* = 7r) ones have nonvanishing expectation value of nematic 
tensor operator Q a p- 



III. NEMATIC ORDER PARAMETER FOR SPIN 
PARTIALLY POLARIZED STATES: 
GENERAL CHARACTERIZATION 

As we have already seen in the previous section, the 
complication of characterizing nematic order when spins 
are partially polarized comes from the explicit symmetry 
breaking induced by external fields. So in this case one 
has to deal with the tensor Q a p which has nontrivial el- 
ements even without nematic order. To distinguish the 
spontaneous symmetry breaking due to the formation of 
nematic order from explicit symmetry breaking in the 
presence of polarization, special care needs to be taken 
of the elements which are induced by spin polarization. 

A general scheme to project out nematic order param- 
eter tensor appears to be possible, similar to what is car- 
ried out in the previous section. Assume spins are po- 
larized along direction s (unit vector). Introducing two 
projection tensors 

1 3 1 

n i/3 = -^ie a /3 7 s 7 , n Q/3 = ^( s « s /3 - g<W), (14) 

we again are able to define a projected nematic order 
parameter as in Eq.(9). 

When the nematic symmetry is broken along the direc- 
tion ui (unit vector), in the large N limit one can easily 
demonstrate that s • ui = following the algebras in Ref . ( 
19, 11); and ui and s always appear to be orthogonal. 
One can further define 



m = s x u>. 



(15) 
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Then oj, m and s form an orthogonal triad. 

It is possible to verify the validity of the definition 
for nematic order parameters given Eqs.(9),(14). For in- 
stance, one can consider the following spin partially po- 
larized nematic states 



l*> = 



[(cos |n + i sin |m) a ^>t ] 

7m 



N 



| vac) 



(16) 



We have assumed that n and m are orthogonal, i.e. 
n - m = 0; £ varies from to tt. Following the discussions 
in Ref.( 19, 11, 13), states specified in Eq.(16) with £ = 
form a complete set of TV-particle-condensate wave func- 
tions. And a condensate with total spin S, S z (S < N) 
can be expressed in terms of spherical harmonics Ys,s z (n) 
in this representation. 

When £ 7^ 0, a state given above is polarized along a 
direction perpendicular to n and m. Indeed, 



(S) = iVsin£n x m, s = n x m. 



(17) 



A direct calculation shows that the projected nematic 
order parameter can be expressed in terms of three or- 
thogonal unit vectors (n, m, s), 



Q 



a/3 



N 



£ £11 

cos 2 -n a n fi + sin 2 -m a m f:} + -s a s - -5 af} . (18) 



Note that the projected nematic order parameter is trace- 
less and fully symmetric 20 . One can further truncate the 
projected matrix in the (n, m)-plane perpendicular to s 
and indeed find that 



n p nm _ ^_ ( cos£ 
Wa P 2 I -cos£ 



(19) 



which is diagonal when a, (3 are chosen to be along the 
axis n or m. 

The projected nematic order parameters Q„p, Qap m 
vanish when spins arc fully polarized or £ = 7r/2; and the 
nematic matrix has zero eigenvalues. When £ ^ tt/2, the 
matrix has nontrivial eigenvalues ±7V/2cos£. The ne- 
matic axis u) therefore lies along the direction of n when 
£ varies from to tt/2 and along the direction of m when 
from tt/2 to tt. The nematic matrix eigenvalues reach 
maxima when £ = or tt, representing spin unpolarized 
nematic states. 

Obviously, nematic symmetry order can develop along 
an arbitrary direction in a plane perpendicular to s. In 
fact, an 0(2) rotation of the orthogonal basis (n, m) 
along s by a <fi angle, while leaves (S) invariant, results 
in a new nematic state with easy axis lu; 



uj = cos <f>n + sin 4>m 
if £ e [0, tt/2], and 

u> = — sin 4>n + cos </>m 

if£e [tt/2, 4 



(20) 
(21) 



In terms of the easy axis u, the projected order param- 
eter can be conveniently expressed as 



Q^ ro = AT|cos£|(u; aW/3 --o a/3 ). 



(22) 



To summarize, we have shown that a projected trace- 
less nematic tensor order parameter should be introduced 
to study nematic ordering in the presence of external 
fields. 



IV. FERROMAGNETIC XXZ MODEL AS THE 
EFFECTIVE HAMILTONIAN CLOSE TO 
CRITICAL POINTS 

A. Phenomenology 

To study the magnetically stabilized nematic order, 
we consider a limit when the exchange interaction J ex is 
much less than E s . For an even number of particles per 
site and in the absence of external fields, the ground state 
is a spin singlet Mott state and nematic order is absent. 
The development of nematic order first occurs when 



2H Z w 6E S > J e 



(23) 



The Hilbert space for the whole lattice is a direct prod- 
uct of spin towers Hk at each site 



Hto = Hi <g> H 2 ® H 3 



(24) 



The on-site Hilbert space Hk is spanned by (N + 1)(N + 
2)/2 states, with spins equal to 0,2,4, N; the dimen- 
sion of the Hilbert space for the whole lattice is 



Vto - ( 



(N- 



l)(N + 2) Vt 



(25) 



where Vt is the number of lattice sites. 

When the external fields satisfying the condition in 
Eq.(23) are applied and when J ex = 0, at each individ- 
ual site the first excited state \S = 2, S z = 2) and ground 
state 15 = 0, S z ) are nearly degenerate and are far away 
from other excited states. At the point when the field 
reaches a value so that 



3E S , 



(26) 



level crossing between the the ground state and first ex- 
cited states occurs as mentioned briefly in the previous 
section. Following discussions in section II, if the hop- 
ping or the exchange energy is set to be precisely zero, 
then magnetization jumps and Q^g = 0. In this case, 
nematic order is not induced by external fields. 

As shown in the previous section, for nematic order to 
be present, two nearly degenerate states have to be in 
a linear superposition. In this sense, naturally it is the 
exchange process in the vicinity of level crossing points 
which eventually leads to nematic order which doesn't 
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exist in zero fields. This observation leads us to truncate 
the on-site Hilbert space into a two-dimensional one for a 
pseudo-spin. The truncated Hilbert space for the whole 
lattice is then a product of pseudo spin Hilbert space Sk 
at each site k 



Hrt = «5i <g> <S 2 ® S 3 ... <g> S k ® ... 



(27) 



and the on-site pseudo spin Hilbert space consists of 
two states: 

IT) = \S = 2,S z = 2),\l) = \S = 0,S x = 0). (28) 

For two particles per site, the microscopic wave function 
of these two states are given in Eq. (7) . 

The dimension of the truncated space Vrt is exponen- 
tially small compared with the original one T> T0 ; i.e., 



V Tt = 2 Vt < £>to; 



(29) 



it is also independent of the number of particles per site. 
The phenomenology for different even numbers of parti- 
cles per site is therefore identical. 

The hopping between two nearest neighbors in lattices 
introduces exchange interactions between pseudo-spins. 
We will present results of microscopic calculations in the 
following subsection. Here we provide a phenomonology 
of this model. To facilitate discussions, we define | f ) and 
I J.) to be two eigen states of Pauli-matrix a z , 



^|T) = IT),^ll) = -U). 



(30) 



Note that these two pseudo spins are also eigen states of 
the spin operator S z . Therefore the pseudo spin algebra 
corresponds to the projection of usual SU(2) spin alge- 
bra in the truncated pseudo spin space. For instance, one 
can verify the following mapping 



h{a z + 1) -» S z , ha + -» S+, tur- 



pi) 



An important and obvious fact is that single parti- 
cle hopping conserves the total spin of two sites and its 
projection along all directions including the z-direction. 
Following Eq.(31), this conservation of spins implies that 
any induced exchange coupling have to as well conserve 
the pseudo spin defined along a z axis in the presence of 
external Zeeman fields. Furthermore, the superexchange 
due to virtual hopping between two bosonic S = 1/2 
pseudo spins results in a ferromagnetic coupling which is 
to be further verified by microscopic calculations. 

Based on the above considerations, one concludes that 
the effective Hamiltonian in the truncated space should 
be 



H, 



eff 



(kl) 



-(/3+l)e £ Pkz&lz — £0 hz^a kz . (32) 

(kl) k 



Here eg, (3 depends on microscopic details of states and 
should be a function of the number of particles per site 
and 7)\, the ratio between E s and E c . h z further de- 
pends on 77 (the ratio between J ex and E s ) and the ratio 
between external fields H z and J ex . 

One can easily recast the Hamiltonian into the follow- 
ing ferromagnetic XX Z model in an effective external 
field along the z-direction, 



H 



xxz 



<Jka<7la 

(kl) 



P^^VkzVlz ~ h z ^ <7 kz 
(kl) k 



(33) 



Because external magnetic fields are applied along the 
z-direction, with the particular choices of eigen states for 
the pseudo spin Pauli matrix a z in Eq.(30), the Hamil- 
tonian in Eq.(33) also has an 0(2) invariance in the xy- 
plane. This 0(2) symmetry represents the 0(2) nematic 
symmetry we are going to examine. The relation between 
the symmetries of the pseudo spin model and the micro- 
scopic model for spin-one bosons has been addressed in 
previous sections. 

In general, the truncation can be applied in the vicini- 
ties of all critical points where level crossings between 
\S, S z = S) and \S + 2, S z = S + 2) occur, S + 2 < N. 
One arrives at the same phenomenology as for the level 
crossing between the first two states. Of course eo, (3 and 
h z then depend on the states involved in level crossings 
and are functions of S, S = 0, 2, 4, ...N — 2. In the next 
few subsections we are going to calculate eo,/? and h z . 



B. Calculations of parameters eo,/3 and h z in the 
XXZ model 

Microscopic calculations of j3 and h z though straight- 
forward are pretty involved. We present results in a few 
limits; detailed calculations can be found in appendix A. 



1. Two particles per site 

There is only one level crossing in this case. One can 
verify that 



i 



(UQM = -if , 



0-7 



. V3 V3 





(34) 



(35) 



and 



(I \Qp a \ T) = (T \Q a p\ l) f , (I \Qpa\i) = 0. (36) 



Using the Hamiltonian in Eq.(2) and taking into ac- 
count of these matrix elements of Q a p in the truncated 
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Hilbert space, one obtains the results for eo, (3 and h z . In 
this particular case, one finds eo = 4/3, (3 = which 
implies an 0(3) symmetry when the effective field h z 
vanishes (but with a finite external Zeeman field H z ). 
It leads to a symmetry higher than the 0(2) one in the 
original problem in the presence of Zeeman field H z . 

We believe that the 0(3) symmetry found for 2- 
particles per site is accidental and can be removed by 
taking into account contributions of order of rji = E s /E c . 
(see Appendix A for details). The final result can be sum- 
marized in the following equation, 



1 + 



h z = -- 



9E„ 



Es 
E c 
3H 



,0 = - 



3E S 



E c + E s 



S T — 2^- T 

°^ ex ^ £ ^ ex 



4J e 



(l + t) 



(37) 



H - E s si - h z s z 

k 

-Jex^2[Quf3(nk)Qp a (ni) + h.c] 

(ki) 



(42) 



where S = in x d/dn, the spin operator is defined as 
the angular momentum of 0(3) quantum rotor. It is a 
conjugate variable to director n, 

[Sa.n^] = -ie Q( 3 7 n 7 . (43) 
Again the matrix elements of Q aj 3 are calculated below 



1 



2 



(t|Qa/j|T>=™ o 2 

- ' 1 -4 



(44) 



2. Four particles per site 

Close to level crossing between |0, 0) and |2, 2), we find 
that 



ii i 

22 

(T 10^10 = ( -i i 







UlO^lt) 



— — 

u u 21 



V 15 'V 15 U 



15 

/ 



V o 



(38) 



(39) 



The corresponding parameters eo, (3 and h z are 

H z 536 



_ 56 351 15 3E S - 

e ° ~ 15'^ ~ ^686' z " ~56 (_ 1Z 



147 



). (40) 



The effective XXZ model has the desired 0(2) symmetry 
in the plane perpendicular to the external field. For four 
particles, level crossing also happens between |2,2) and 
|4,4) states. Similar calculations have been carried out 
and presented in Appendix A. 



aio^it) 




(45) 



The matrix element Q xy vanishes in Eq.(44) as an ar- 
tifact of the large- TV approximation. One then obtains 
all parameters for the XXZ effective model, 



2 39 15.3£ s - 

e °-T5' /3 = ^49'^ = "Y (_ 1Z 



147)- ( 46 ^ 



It is possible to generalize this analysis to level cross- 
ing between high spin states \S, S) and \S + 2, S + 2) 
(S < N — 2) and results are qualitative the same (see 
Appendix B). In all cases, eo is positive and is nega- 
tive. In the next section, we are going to examine the 
consequency of this model. Particularly we investigate 
the implications on magnetically stabilized nematic or- 
der and physics around critical points. 



V. PHASE BOUNDARIES OF THE XXZ MODEL 
AND HOLSTEIN-PRIMAKOV BOSONS 



A. Phases of XXZ model 



3. Large-N limit (even N) 

In the large- N limit, one can describe the collective 
ground state and excited states in terms of spherical har- 
monics in a quantum rotor representation. So the spin 
singlet ground state and polarized \S = 2, S z = 2) wave 
functions are 

IT) = i/fsin^ex P (^), 1 1) = J=. (41) 

In the quantum rotor representation, the Hamiltonian 

:„19,10,11,13 



The general phase diagram in the (j3, ft. z )-plane can be 
easily obtained first in a mean field approximation. Later 
on we argue that the phase boundaries and solutions ob- 
tained in this way in some part of the plane are actually 
exact (see Fig. 2). In the mean field approximation, 
we introduce s as a unit vector order parameter which 
defines the orientation of spin, 

(a) = 2Ss; (47) 

here S = 1/2 is the pseudo spin. 

The s-dependence of the total energy comes entirely 
from the terms proportional to (3 or h z ; that is 
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E 



= const - 4df3S 2 s 2 z - 2h z Ss z . 



toJexVr 



(48) 



where d = 3 is the dimension of three dimensional cubic 
lattices. s 2 varies from —1 and 1. Minimizing the en- 
ergy with respect to s one obtains mean field solutions 
for various ground states. 

Following Eq.(48), when h z > -2d(3 and h z > 0, the 
mean field solution is 



s= (0,0,1) 



(49) 



representing a fully polarized state which we call Up- 
Polarized (UP) phase. 

When h z < 2d[3 and h z < 0, the mean field solution is 



s= (0,0,-1) 



(50) 



representing another fully polarized state which we call 
Down-Polarized (DP) phases. 

In addition, when — 2d(3 > h z > 2d(3 and (3 < 0, the 
mean field solution is 

s = (sin 6 cos $, sin 9 sin $, cos 9), cos 9 = (51) 

representing a ferromagnetically ordered (F0) phase 
which breaks the in-plane 0{2) symmetry spontaneously. 
Solutions of different angle <& are degenerate and the vac- 
uum manifold is a unit circle S 1 . 9 varies from it to 
when h z increases from 2d(3 to —2d(3. 

Three phases are separated by a first order transition 
line along the h z — axis which starts at point (0, 0) 
and ends at (oo,0), and two other second order phase 
transition lines (see Fig. 2). These two lines arc defined 

by 



2d(3±h z = 0; 



(52) 



both terminate at point (0,0). Finally (0,0) is a tri- 
critical point. 

Along the first order transition line, the UP and DP 
states become degenerate and the ground state breaks 
Z2 or Ising type of symmetry spontaneously. At the 
tricritical point (0,0), the XXZ model is 0(3) rotation 
invariant and the ground state breaks 0(3) symmetry 
spontaneously At this point, UP, DP and F0 states are 
all degenerate. 

In UP and DP phases, the microscopic wave functions 
for ground states are, respectively, 



iffT>=nit>*.i0i>=n u >*; 

k k 

Cfcz|T)fc = \\)k,^kz\i)k = -|l)fe- 



(53) 



In the 0(2) ferromagnetic phase, 
\9f) = J||0) fe , 

k 

|0) fe = cos ® cxpH$)|T) fc + sin ®||) fe ; (54) 



and 9 is a function of (3,h z , cos9 = ^fe, \h z \ < 2d\/3\. 
Solutions in Eq.(54) are degenerate in the ^-manifold 
where exp(z$) lives and represent spontaneous 0(2)- 
symmetry breaking states. 

By examing the microscopic wave functions of UP and 
DP states, we notice that the UP and DP states are non- 
degenerate exact eigenstates of the pseudo spin operator 
S z = J2 k <Jk z - Meanwhile, the total pseudo spin pro- 
jected along the z-axis is a conserved quantum number. 
So these UP and DP solutions are exact eigenstates of the 
XXZ-Hamiltonian. In the next subsection we are going 
to show that they are actually exact ground states when 
(3 is positive; furthermore we argue that UP or DP states 
are also exact ground states even when (3 is negative and 
h z > -2d/3 or h z < 2d(3. 



B. UP and DP states as exact ground states 



When (3, h z are both positive, the UP state presented 
above is a ground state of both the 0(3) isotropic term 
in the ferromagnetic XXZ model and the terms involv- 
ing (3, h z . So naturally the UP state is the exact ground 
state of the XXZ model in this limit. Similarly when f3 
is positive and h z is negative, the DP state is the exact 
ground state. 

When (3 < but outside the triangular defined by the 
two critical lines 2d(3 ± h z = 0, we are not able to prove 
rigorously that eigenstates in Eq.(54) are exact ground 
states. However, we would like to show that they are 
locally stable and therefore we argue that they are likely 
to be the exact ground states. 

To carry out this part of discussions, we study the 
magnon excitation spectrum in UP and DP phases and 
show that one-particle magnon excitations are also ex- 
act eigen states; furthermore they are fully gapped ex- 
cept along the second order transition lines. The most 
straightforward approach to study these excitations is to 
use the Holstein-Primakov boson representation for the 
XXZ model. 

In the Holstein-Primakov representation, all spin op- 
erators are expressed in terms of Holstein-Primakov 
bosons, 



2S 



C'C) c 



a 



2(S-c t c) 



(55) 

(56) 
(57) 



(c) is the creation (annihilation) operator of bosons sat- 
isfying the usual bosonic commutation relations [c, c 1 ] = 
1; and the raising and lowering operators are defined as 



(58) 
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One can furthermore verify that 

WatVp] = *2e a /3 7 cr 7 , er • a = 4S(S + 1) (59) 

The Hamiltonian of the XXZ model then transforms 
into: 



H 



XXZ 



-2 V t 25 - c\ c k\ck4^2S - c\ Cl 
(kl) 



(u) 

(kl) 

+2(h z + 45(1 + p)d)J2c{c k . 

k 



(60) 



In deriving Eq.(60), we have neglected a constant term 
— 4(1+(3)S 2 — 2h z S for each lattice site. In a semiclassical 
approximation, one indeed recovers the results obtained 
in the previous section. Again S = 1/2. 



(-1,0) 
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ii 



FIG. 2. Phases in the ferromagnetic XXZ model. /3 
varies from —1 to +oo. Region 1,11 and III represent 
Up-polarized(UP), Down-Polarized (DP) and Ferromagnetic 
Ordered (FO) phases respectively. Along the blue lines 
(—2d/3 ± h z — 0), transitions are continuous while along the 
red line the transition (h z = 0, j3 > 0) is a first-order one. 
Point (0,0) is the 0(3) symmetric tricritical point of the fer- 
romagnetic XXZ model. As q goes to zero, the interactions 
between magnons are repulsive when (3 < and attractive 
when f3 > 0; along the dash line (/3 = 0), magnons are 
non-interacting. The solutions in the shaded region can be 
obtained in a dilute gas approximation. 

To study the excitation spectrum in region I (or II), we 
first examine the Hamiltonian in Eq.(60) in a one-particle 
subspace next to the exact eigen states of UP (or DP). 
First one notices that a UP state is an exact vacuum for 
Holstein-Primakov bosons; that is 



Cfclffr) = °>vk Z \gr) = \gt)- 

for any lattice k. 



(61) 



One-particle excitations we are interested in live in a 
subspace of single Holstein-Primakov boson; that is in a 
space spanned by states 



4l5T> 



(62) 



defined at each lattice site k. Since the total number op- 
erator of Holstein-Primakov bosons commutes with the 
Hamiltonian 



[JV C , H] = 0, N c 



J2 c i Ck > 

k 



(63) 



N c is a conserved quantum number. We can then di- 
agonalize the Hamiltonian in this one-particle subspace 
where iV c = 1. 

In the subspace, we obtain the following effective 
Hamiltonian 



rro.p. 

n xxz _ 
e q = 8S(d 



E 



e qCqC q ; 



(64) 



The superscript o.p. stands for the "one particle" sub- 
space. Eq.(64) indicates the dispersion relation of one- 
particle states. 

In particular, it yields a fully gapped magnon spectrum 
in region I; the gap vanishes only along the second order 
transition line where (3 < and h z = ±2df3. Especially 
magnons are fully gapped along the first order phase 
transition line j3 > and h z = 0. When \h z \ < —2df3 and 
[3 < 0, one-particle states (or magnon excitations) start 
to have lower energies than the vacuum state. This in- 
dicates condensation of Holstein-Primakov bosons which 
we are turning to. 

To conclude we find that DP and UP states are ex- 
act ground states of the XXZ model in region I and 
II; magnon excitations in these phases are fully gapped. 
Along the mean field second order transition lines, 
magnons become gapless excitations. Further decreas- 
ing h z results in instability of magnon excitations. So we 
believe that the transition lines in the mean field theory 
represent the exact phase boundaries. 



C. Condensation of Interacting Magnons and 
Emergence of Ferromagnetic Order in the XXZ 
model 

As discussed in the introduction, the dynamics of con- 
densation of magnons depends crucially on the interac- 
tions between magnons. To study the region close to 
critical lines where the condensed particle density should 
actually be very low, we only take into account two-body 
interactions and apply a dilute gas expansion. The re- 
sults we derive in this subsection are valid in the shaded 
critical regions (see Fig. 2) where 







\h z ± 2d/3\ 



< 1. 



(65) 



And as /i z ± 2d(3 approaches zero, the results become 
exact 21 ! We will present the calculations along the up- 
per transition line defined by h z + 2d(3 = 0; the results are 
then generalized to the lower transition line h z — 2d(3 = 0. 

In the dilute gas limit which interests us, the number 
of Holstein-Primakov bosons per lattice site is much less 
than one i.e. 



n c = (c{c k ) < 1. 



(66) 



For this reason, one can expand the nonlinear operators 
of a in terms of n c ; especially, 

V2S^=V2s(l-^ + 0(n 2 c )..)j (67) 

for S = 1/2. The resultant many-body Hamiltonian up 
to the second order of n r is 



H 



xxz 



— J]] e q c q c q 



-Ad/3 ^ ci,_ l _ n „ct,„_„,C a ,<a 



qi +q3 q2 - q;? °qi ^qs • 



(68) 



qi,q2,q3 



The first term is identical to H°xxz-> * ne exac t Hamil- 
tonian projected in the one-particle subspace; the sec- 
ond term describes magnon-magnon interactions. This 
Hamiltonian is applicable in a dilute limit where Eq.(66) 
is satisfied. The sign of interaction at q — or small 
q limit is determined by (3. When (3 is positive, magnon 
interactions are attractive and when negative magnon in- 
teractions are positive. 

When magnons are ideal {(3 = 0), all magnons con- 
dense when the energy gap in the spectrum closes at 
h z = 0. This leads to an abrupt jump in magnetization 
which corresponds to the field-driven first order phase 
transition from UP to DP phase at the tricritical point 
(0,0) (along the dashed line). And in this case external 
fields do not induce nematic order. This is consistent 
with mean field results discussed in the previous subsec- 
tion. One can in principle generalize this argument to 
the case when (3 < and arrive at similar conclusions. 

When magnons' interactions are positive, following 
Eq.(68) the chemical potential of magnons in the dilute 
gas limit (differing from no OI atoms) is 



// 



-8n d(3 + O(n 2 ) 



(69) 



where no is the number of magnons per lattice site. This 
is similar to weakly interacting gases of bosons in contin- 
uum limit 22 . The chemical potential defined in this way 
only depends on intrinsic parameters (3 which have been 
evaluated microscopically and is independent of external 
Zeeman fields. The energy of the magnon condensate per 
lattice site is therefore 



E(n ) 



2{\h z \ +2df3)n - 4nld(3. 



Minimizing the energy with respect to no yields 



n 



2 [ 2<2/3 J 



(70) 



(71) 



which is a continuous function of h z . n is much less 
than one in the critical region where Eq.(66) is satisfied. 
At the transition point \h z \ = —2d(3, the magnon den- 
sity per lattice site either vanishes or is equal to one, 
i.e. no = 0, 1. Furthermore, if one extrapolates to the 
h z = case, one obtains no = 1/2, that is half magnon 
per lattice site. a z — as expected. 

Note that the ground state in this case is not the vac- 
uum of Holstein-Primakov bosons but instead the vac- 
uum defined by Bogolubov quasi-particles. The Bogol- 
ubov excitations are created by the following operators 



7 q = "(q)c q + v{q)c 

« 2 (q) = J(i 



+ n 



v 2 (ci)=--(l- 



^/e 2 +2e^u 

e q + M \ 



y / e q + 2e q /x 



(72) 



where u is the chemical potential of magnons defined be- 
fore and the kinetic energy e q = 2|q| 2 a 2 is written in 
a dimensionlcss unit. The dispersion of quasi-particles 
is phonon-like uo q — v s \q\ at small energies; taking into 
account the chemical potential in Eq.(69), we obtain 



v s = v s0 1 



\h z 



1 + 2c[j3 ,Vs0 = ^V 3 -^- 



(73) 



This agrees with the semiclassical solutions obtained in 
Eq.(C12) in Appendix C. 

These results indicate that the physics in the XXZ 
model close to the second order critical lines is in- 
deed equivalent to interacting dilute magnons defined 
by Holstein-Primakov bosons; especially the emergence 
of ferromagnetic ordering occurs when condensation of 
magnons takes place. As in the usual c-number approxi- 
mation for condensed bosons, we approximate 



c Lo = c q=o = V n oV T exp(i^). 



(74) 



Substituting this result into the expressions for g ± , a z in 
Eq.(57), we obtain 



(75) 



(a x ) = 2v2Sy / n ~ cos $, 
(<T y ) = 2\/25 v / no"sin$, 
(a z ) =2S-2n . 



And again 5" = 1/2. 
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Correspondingly, the Bogolubov quasi-particles repre- 
sent the spin wave excitations in 0(2) ferromagnets. Fol- 
lowing Eqs. 57, 72, one can express a XtV in term of 7^ 
and 7, 



5<r x {r) = ^2 cxp(iq • r) 

q#0 



2y/V T 
2S 



(u(q) - u(q))(7q + 7-q)> 



5(J v( V ) = J2 eXP (* q ' ^ O- /T^ ("( q ) + U (q))(7q _ 7-q)- 
q^O 2l ^ VT 

(76) 

Consider a single quasi-particle state 

|q )=7q l^c). (77) 
Spin correlations in this single particle state are 



2S 



((6a x (r) - 6a x {0)) 2 ) = —(u(q ) - u(q ))^sin 



s 2 ■ 2 qo • r 



V T 
2S 
Vt 



((Sa y (r) - 5a y (0)Y) = ^(u(q ) + «(q )) 2 sin 2 S£_£, 



S 



((Sa x (r)Sa y (0))) = -— sinq • r. 



(78) 



Remarkably, the corresponding orientation of pseudo 
spin s derived in the dilute gas approximation is pre- 
cisely the same as the semi-classical results obtained in 
subsection A; close to the critical line, we notice that 



cos = 1 — 2n 0} sin 6 = 2-y/no. 



(79) 



In the next section we are going to discuss the impli- 
cations of the mapping on correlated atoms, especially 
magnetically stabilized nematic order. Since the semi- 
classical solutions turn out to be exact along the critical 
lines, we would like to believe that solutions are also valid 
in the ferromagnetic ordered phase, at least qualitatively 
and can be extrapolated deep into that phase. 



VI. NEMATIC ORDER AND PHASE 
BOUNDARIES OF MAGNETICALLY 
STABILIZED NEMATIC MOTT STATES 



Let us turn to the problem of Mott states of spin- 
one bosons. Following discussions in section II, one finds 
that DP states correspond to spin singlet Mott (SSMI) 
states and UP states to spin fully polarized Mott (SF- 
PMI) states. 

The FO states breaking the 0(2) symmetry represent 
quantum spin nematic states with easy axis determined 
by the projection of pseudo spin orientation s in the xy- 
plane. Indeed, for two particles per site the wave function 
of FO states in Eq. (54) indicates the following spin cor- 
related Mott states for spin one bosons, 



NM 



2V2 



■ sin — exp^ — ) Hot Lr c 
2 ^ 2' V6 



}\vac). 



(80) 



is a function of [3 and h z as given in section V.A, 



e 1 h z . e 1 h z 

cos — = \ , sm — = \ — ; 81 

2 y 2 Adj3' 2 y 2 4d/3' V ' 

As h z varies from —2d(3 to 2d(3, varies from to tt. 
And $ G [0,27r] represents an ^-manifold for the spon- 
taneous symmetry breaking solutions. 

For cold atoms, our calculations show that in all cases 
/3 is negative. When magnetic fields are varied, the tra- 
jectory in the (3 — h z planes (see Fig. 2.) moves vertically 
at a given negative j3. And as magnons are repulsive, 
magnetic fields stabilize nematic order via the continuous 
process of condensation of magnons. So as the magnetic 
field increases, phases encountered are spin singlet Mott 
states, nematic Mott states (partially polarized) and spin 
fully polarized states. Here we will focus on the nematic 
state. 



Sin 




FIG. 3. Magnetization Af z (in units of h) and sinO as a 
function of h = h z /2d(3. sinO defined in Eqs. (51), (81) is pro- 
portional to the nematic order. Results around upper or lower 
critical points are obtained in a dilute gas approximation and 
are exact. Notice that the nematic order sinO reaches the 
maximum at zero h z where M z is equal to h. 

The projected nematic order parameter for the con- 
structed nematic Mott state is given in Eq.(ll) with 
6 = $ and 6 = 0: 




h 2 1 



w = (cos-, sin-). 



(82) 



The nematic order vanishes along the second order crit- 
ical lines \h g \ = —2d(3 and reaches maxima when level 
crossing takes place in an isolated lattice site, i.e. at 
K=0. 
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This spin partially polarized nematic Mott state (SPP- 
NMI) has spin polarization 



(83) 



Spins are fully polarized at one of the critical lines 
(M z = 2h) when h z + 2d(3 — and the spin polarization 
vanishes (M z = 0) at the other critical line h z — 2d(3 = 0. 
In between, M z varies continuously from to 2Ti and is 
precisely equal to h when h z vanishes and the nematic 
order reaches the maximum. 

Finally, the phase boundaries for SSMI, SFPMI and 
SPPNMI can be obtained by substituting the field depen- 
dence of (3, h z derived in section IV. B into the expression 
for critical lines in the XXZ model, 



hz (^,E±)±2dP& = 0. 

'J ex J ex J^c 



(84) 



This results in critical fields for various N. Especially 
one determines the upper and lower critical fields (H^ c ) 
between which nematic order develops for the first time 
when magnetic fields increase from zero. 

For N = 2, the upper and lower critical fields are 



8 E 
H zc — 3E S — -J ex ± 24— J ex ; 



for N = 4, these fields are 



Ht c = 3E s -™J ex ± 2 -^lj ex . 

zc s ex ex 



(85) 



(86) 



At the large- N limit, one obtains 

H zc = 3E S - —J ex ± —J ex . (87) 
We have set d = 3 in deriving Eqs.(85),(86),(87). 



VII. EFFECTS OF QUADRATIC ZEEMAN 
COUPLING 

In this section, we are going to briefly discuss the effect 
of quadratic Zeeman terms which generally are present 
in atomic gases 16 ' 17 . This kind of external perturbations 
only conserves the spin projection along the direction of 
external fields but does not conserve the total spin of 
the many-body states under consideration and therefore 
has distinctly different effects on spin singlet Mott states. 
Namely, such external fields would induce nematic order 
at any small but finite coupling. In other words, spin 
singlet Mott states are unstable with respect to these 
perturbations. 

To demonstrate this phenomenon, we consider spin sin- 
glet Mott states in the presence of the following quadratic 
Zeeman perturbation 15,19 , 



Hp = -H q . z . ^2 Qa(i{k){n a n ( j - ^5 a p); 



(88) 



here Q a p(k) is the nematic operator defined at the be- 
ginning of section II (Eq.2) and n characterizes the ori- 
entation of quadratic Zeeman fields which we choose to 
be along the z-direction. H q . z . is the strength of the 
quadratic Zeeman coupling. (This perturbation differs 
slightly from the one used in Ref.( 15, 19) by a singlet 
operator which doesn't contribute to the quantity we are 
calculating here.) 

One notices that indeed this quadratic Zeeman term 
does not communte with the total spin operator defined 
at any individual lattice site, Sf; however it does com- 
munte with the operator Skz, 

{^kzi'^Qa' fi'{k')S a > z 5 (jl z ] = 0, 

k' 

Skz = -i^e za 0tpl a tp k f3- (89) 
k 

So what it does is to cause transitions between states 
with different on-site spin quantum numbers Sk but with 
identical spin projection along the z-direction, Skz', thus 
it does not lead to transitions between different Sfc z sub- 
spaces. For a spin singlet Mott state, this perturba- 
tion results in transitions between on-site singlet states 
(Sk = 0) and non-singlet states (Sk ^ 0) in the subspace 
of S fez = 0. 

For two particles per site and in the zero hopping limit, 
we find that these transitions lead to coherent superpo- 
sition of two states in S z = subspace: |0,0) state and 

|2,0) = ^(3VM - VM)Nc),a = x,y,z. (90) 

For instance, in the first order perturbation expansion 
the ground state wave function is 



(5* 



V2H b . z . 
9E S 



(91) 



A direct calculation shows that the nematic order is in- 
duced continuously as the quadratic coupling is applied 



(92) 



This dependence is very different from the linear-Zeeman 
field dependence of nematic order which exhibits a criti- 
cal field below which spin singlet Mott states are stable. 
As expected, quadratic Zeeman effects are more effective 
in stabilizing spin nematic Mott states. Furthermore, the 
easy axis of the induced nematic tensor order parameter 
is pinned along the direction of external fields, n and the 
resultant states are Ising nematically ordered instead of 
0(2) or 0(3) nematic states discussed before. 
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VIII. CONCLUSIONS 

To summarize, in this article we have investigated mag- 
netically stabilized fluctuating spin nematic order. We 
have shown that nematic order can develop when two 
non-nematic states at a lattice site are in a linear super- 
position in the presence of external fields. When exter- 
nal fields are applied, even small superexchange coupling 
could lead to such a linear superposition and nematic 
order emerges even though no spontaneous symmetry 
breaking occurs in zero field. 

We have also mapped the problem of spin-one bosons 
with antiferromagnetic interactions in an external field 
to the ferromagnetic XXZ spin (S = 1/2) model. We 
find that the field-driven quantum phase transitions be- 
long to the universality class of the ferromagnetic XXZ 
model (S = 1/2). Spontaneous symmetry breaking in 
the xy-plane in this effective ferromagnetic XXZ model 
corresponds to planar nematic ordering in the underly- 
ing atomic states. In all non-nematic Mott states which 
interest us, interactions between magnons are repulsive. 
Therefore when the external field reaches a critical one, 
condensation and thus phase transitions are continuous. 

We also show that the breaking of the U(l) symme- 
try in magnon Bose condensates results in breaking of 
the 0(2) nematic symmetry in the xy-plane perpendic- 
ular to external fields. The Bogoliubov quasi-particles 
of condensates are precisely the spin wave excitations in 
the 0(2) nematic states. So the nematic order is stabi- 
lized when the field exceeds a critical one and magnons 
condense. We have also obtained the microscopic wave 
functions of ordered states and spin wave excitations. 

Finally we find that for a spin singlet Mott state the 
fluctuating nematic order can be stabilized by any small 
but finite quadratic Zeeman effects. Namely, the nematic 
order parameter varies continuously in the presence of 
quadratic Zeeman effects. 
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APPENDIX A: EFFECTIVE XXZ MODEL FOR 
VARIOUS NUMBERS OF PARTICLES PER SITE 

1. Microscopic Hamiltonian 

To break the 0(3) symmetry in the XXZ model for 
two-particles, we keep higher order terms of o(E s /E c ). 
The effective Hamiltonian for spin-one bosons in the 



Mott state of an optical lattice in the presence of a mag- 
netic field (in the z-direction) can be derived as: 

H=(e s - 6J ex ^j ]T S 2 fe -H z J2k (Al) 



-2J ex 1 + 



El 



a 

(ki) 



E : 



+2J e ^^(S fe + S J ) 2 



(ki) 



E 

+j ex ^J2 T ^ AkBi + AiBk ^ 

c - (hi) 



(A2) 



We have introduced the operators B v ^ = ip^ip^ and 



2. N = 2 case 



For N = 2we have only the possibility of looking at the 
transition between the states | j) = |2, 2) and | J.) = |0, 0) 
. The relevant non vanishing matrix-elements are (again 

(TIQa/jU) = (TIQ ct/3 |T> t ): 



(TIQa/sIT) 
UI<Mt) 



1 

3 


i 







—i 


i 

3 





)■ 








2 
3 




i 




i 










I 






:) 


V3 




V3 














1 -i 

(1\A\l) = V3 \ -i -1 



-1 -i 

(l\A\1) = V3 \ -i 1 





(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
(A8) 

(A9) 
(A10) 
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Using these results we find that the effective Hamilto- 
nian turns out to be: 



3 V E c 
E 



(All) 



(ki) 



Hz - 3E S + -J ex + -—^-Jex ) ^2 a k- 



3 3E, 

This is the Hamiltonian for the XXZ-model: 
Hxxz 

with 



■5> fc -<r,-/?5>fo*-/i,5>2 (A12) 



(ki) 



(ki) 



eo = 3 + 



0-- 



3E S 



E c + E s 



QE S 3>E[ Z SJex ^ J ex 



4J B: 



(' + *) 



(A13) 
(A14) 

(A15) 



b. |2,2> — |4,4> 

We define: | |) = |4,4) and 1 1) = |2,2). The relevant 
non vanishing matrix-elements arc: 



(TIO^IT) 



UIQ«/»lt) = 



UIQ^U) = 



| 2i 
-2i | 



0- 

/ A /IJ ^ o \ 





V o 



-i A 
1 21 

0- 



o oy 
o 





22 



(A21) 



(A22) 



(A23) 



Using this we get again an effective Hamiltonian in the 
form of an XXZ-model. The parameters are: 



P 



g= ^2MJex = _n 
49 J 84 



h z 



1 



7E S - H z 



344 
49" 1 



(A24) 
(A25) 
(A26) 



3. N = 4 case 



For four particles per site, there are two possible tran- 
sitions: |0, 0) — > |2,2) — > |4,4). We will consider both 
transitions. In both cases the correction of order J ex %- 
turns out to be not particularly interesting. 



a. |0,0) - |2,2) 

We define again | |) = |2,2) and | J.) = |0,0). The 
relevant non vanishing matrix-elements are: 



[\\(ht\\) = \ -\ 4 





22 
" 21 



V 15 V 15 



aiQ^it) 



15* V 15 U 



V o 



14 

o 0/ 



(A16) 



(A17) 



This gives rise to an XXZ model with the following pa- 
rameters: 



56 

60 -15' 



(3 = 
h z = 



351 
686' 
1 

£()Jex 



3F - H 53 S 
0-C/ s n z \/^J x 



(A18) 
(A19) 
(A20) 



4. iV = 3 



For three particles per site we have the transition be- 
tween |t) = 1 3, 3) and | J.) = |1,1). The non- vanishing 
matrix-elements are: 



(T|Qa/j|T> 




(A27) 



(A28) 



(A29) 



The parameters of the XXZ-model in this case are: 
12 



e ° = T 

h z = — 5E S - H z — J ex 



(A30) 
(A31) 
(A32) 
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APPENDIX B: EFFECTIVE HAMILTONIAN FOR 
LARGE N 



APPENDIX C: THE HOLSTEIN-PRIMAKOV 
BOSONS REPRESENTATION 



For large N, N even, we can study all transitions from 
|5, 5} -> |5 + 2, 5 + 2). These states are given by: 



\S,S) = \Y S3 {n)) = {-l)sJ 2S + 1 



4tt 2 2S (5!) 2 



iS +sm s 



(Bl) 



1. Holstein-Primakov bosons in UP or DP phases 

The Hamiltonian of the XXZ-model is given as 



H 



xxz 



The hopping term in the Hamiltonian is just equal to: 
-2J ex ^(n fe -n ; ) 2 (B2) 

(H) 

Introducing | ]) = \S + 2, S + 2) and 1 J.) = |5, S), we 
get the following non-vanishing matrix-elements: 



^ &ka ■ Via ~ P 51 Okz ■ <Jl z - h z ^ G kz . (CI) 

(hi) (kl) k 



(t | fc <T Kn fe .n z ) 2 | T)fe| T>i 

(T |fe(l |((n fe • n z ) 2 | t)fe| I) i 

(T UU | ; (n fc - ni ) 2 | |)fe| T)i 

(I |fe(l |;(n fe • n z ) 2 | J.)fe| |); 



19+ 125 + 25 2 
(7 + 2S*) 2 
7 + 85 + 25 2 
(3 + 25)(7 + 25) 
(2 + 5)(l + 5) 
(3 + 25)(5 + 25) 
3 + 45 + 25 2 



In this subsection, we are interested in Region I (See 
Fig. 2) where h z + 2dj3 > 0. After Fourier transforming 
and in terms of H.P. bosons the Hamiltonian can written 
as 

Hxxz = # (0) + H (2) + # (4) + G((c (t) ) 6 ). (C2) 
Here 

# (0) =~V T (d(J3 + 1) + h z ); (C3) 



(3 + 25) 2 

The effective Hamiltonian turns out to be 



4(1 + 0)d + 2h z - 4 53 cos (la a ) 



CqCqJ 



(C4) 



H, 



eff. 



<feo 



(3 + 25)(5 + 25) 
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<fe;> 



(3 + 25) 2 (7 + 25) : 



^k a i 



and the fourth order term is 



Jex^O Vt 



Vrr, ^ H 



c q 2 c q3 c qi+q2-q3 



qiq2q3 



- (25 + 3)E S -B z - 



8(3 + 5)(l + 25) 
(3 + 25) 2 (7 + 25) 2 ' 

This is clearly an XXZ model with: 
(2 + 5)(l + 5) 



J ex ) 5^°A 



eo = 



(3 + 25)(5 + 25) 



-1 + 



12(5 + 25) 



((7 + 25) 2 (6 + 135 + 95 2 + 25 3 ) 



(2S + 3)E S -B Z 



8(3 + 5)(l + 25) 
(3 + 25) 2 (7 + 25) 2 



(B3) 
(B4) 

J PA 



It is clear that (3 is negative for all 5 and it get's more 
negative if we increase 5. 

For the lowest transition: |0, 0) — > |2, 2) this gives the 
numbers: 



^[exp(-iq 2 „a) + cxp(-zq 3Q a) + 
exp(-i(-qi q 2 + q 3 ) Q a) + cxp(zq lQ a)] 

-4(1 + /?) 53 

c <n c <i2 c qs c qi +q2 -q3 

qiq2q3 

x 53 cx p(*(qi - q3) Q a)- 



Here a is the lattice constant. 

Following Eq.(C4), the energy of the quasi-particles is 
given by 



e q = 4(1 +(3)d-4 53 cos(q a a) + 2h z 



(C5) 



eo 



15 



39 
49 



h z = - 



3E S — H z 



147 



J e; 



(B5) 
(B6) 
(B7) 



a— x,y,z 

where the energy gap in the spectrum is given as 

A(0,h z )=4(3d + 2h z . (C6) 

The fourth order term describes interactions between 
magnons. Indeed, in the small |q| limit the hamiltonian 
can approximately be written as (up to a constant) 



15 



H 



xxz 



£ q C qCq 4/3rf C qi C q 2 C q3 C qi+q: 



2-q3 • 



(C7) 



eq.(C12) agrees with the results derived in the dilute gas 
approximation in section V.C; close to critical lines, one 
notices that sin = and d{3 sin 2 6 = 4d(3no- 



When (3 > interactions between the magnons are at- 
tractive and when j3 < interactions are repulsive. 

To derive these results, we have used the dilute gas 
expansion expansion 

V^Jl= (l-^-I( C t c) 2 + ...). (cg) 



2. Holstein-Primakov Bosons in Ferromagnetically 
Ordered Phase 



The most convenient way to study HP bosons in region 
III is to introduce the following rotation: 



x \ ( cosQx' + sin Qz' 

y) = \ y' 

z } \ cosQx' — sinOz' 



(C9) 



In the semiclassical approximation, by minimizing the 
energy with respect to 6, one obtains the ground state 
solution with cos 6 = — 7^5. 

Consider an expansion over this solution. We get the 
following lowest order terms: 



H 



(2) 

XXZ 
Co J ex 



= ^ (Ad + Adf3 cos 2 9 + 2h z cos 6) c q c q (CIO) 
q 

— ^ (4 + 2/3 sin 2 0) ^ cos q Q a c q c q 



(3 sin 2 <d cos q Q a (c q c q + c q c, 



When 6 = 0, one recovers the results in section V.B. 

Taking into account cos Q = — -J^ in the ferromagnetic 
phase, in the long wave length limit one further simplics 
the result to 



H 



(2) 

XXZ 



4(d — cos q Q a)c q c q 



(Cll) 



-^2 (2d/3sin 2 6) c q c q 
q 

- sin 2 9 (c q c q + c q c q ) . 

q 

This yields the following dispersion 



w q = 2v / 2av / -d/3Wl 



/i 2 



(C12) 
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